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Abstract. CIT 3 is an oxygen-rich long-period variable evolving along the Asymptotic Giant Branch and is one of the most 
extreme infrared AGB objects. Due to substantial mass loss it is surrounded by an optically thick dust shell which absorbs 
almost all visible light radiated by the star and finally re-emits it in the infrared regime. We present the first near infrared 
bispectrum speckle-interferometry observations of CIT 3 in the J-, H-, and A''-band. The J-, H-, and if -band resolution is 
48 mas, 56 mas, and 73 mas, resp. The interferograms were obtained with the Russian 6 m telescope at the Special Astrophysical 
Observatory. While CIT 3 appears almost spherically symmetric in the H- and /("'-band it is clearly elongated in the J-band 
along a symmetry axis of position angle —28°. Two structures can be identified: a compact elliptical core and a fainter north- 
western fan-like structure. The eccentricity of the elliptical core, given by the ratio of minor to major axis, is approximately 
£=123 mas/154mas=0.8. The full opening angle of the fan amounts to approximately 40°. Extensive radiative transfer calcu- 
lations have been carried out and confronted with the observations taking into account the spectral energy distribution ranging 
from 1 fim to 1 mm, our near-infrared visibility functions at 1 .24 ^m, 1 .65 and 2.12 /im, as well as 1 1 /xm ISI interferometry. 
The best model found to match the observations refers to a cool central star with Tcff — 2250 K which is surrounded by an 
optically thick dust shell with r (0.55/im) = 30. The models give a central-star diameter of O, — 10.9 mas and an inner dust 
shell diameter of Oi — 71.9 mas being in line with lunar occultation observations. The inner rim of the dust-shell is located 
at ri = 6.6-R, and has a temperature of Ti = 900 K. The grain sizes were found to comply with a grain-size distribution 
according to Mathis et al. (1977) with n(a) ~ a~^'^, and 0.005 /im < a < 0.25 /im. Uniform outflow models, i.e. density 
distributions with p ~ turned out to underestimate the flux beyond 20 /im. A two-component model existing of an inner 

uniform-outflow shell region (p ~ l/r'^) and an outer region where the density declines more shallow as p ~ l/r'""' proved to 
remove this flux deficiency and to give the best overall match of the observations. The transition between both density distribu- 
tions is at r2 = 20.5ri = 135.77?, where the dust-shell temperature has dropped to T2 = 163 K. Provided the outflow velocity 
kept constant, the more shallow density distribution in the outer shell indicates that mass-loss has decreased with time in the past 
of CIT 3. Adopting Vcxp = 20 km/s, the termination of that mass-loss decrease and the begin of the uniform-outflow phase took 
place 87 yr ago. The present-day mass-loss rate can be determined to be M — (1.3 — 2.1) ■ 10^^ Mo/yr for d — 500 — 800 pc. 
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AGB and post- AGB — Infrared: stars 



1. Introduction 

Within the framework of the Two Micron Sky Sur vey o r 
Infrared Catalogue (IRC, Neugebauer & Leighton 1969), 
UMch et al. ( |1966| ) pubHshed a Hst of 14 bright (at 2.2 fim) 
and very red objects being of potential interest for the 
study of cool and dust-enshrouded stars. Among these so- 
called CIT objects is CIT 3 (= WXPsc = IRCh-10011 
= IRAS 01037H- 1219), an oxygen-rich long -period variable 
evolving along the Asymptotic Giant Branch (AGB) and one 
of the most extreme infrared AGB objects. 



CIT 3 belongs to the very late-t ype A GB stars ha ving 



spectral type of M9-10 (Dyck et al. |1974 Lockwood |1985[ ). 
It is very faint in the optical {V > 20 mag) but bright in the in- 
frared (K=2mag). Given this spectral type and the exception- 
ally large color index ofV — K <; 18, its effective temperature 
can be estimated to be ~ 2500 K or l ess consider ing the effec- 
tive temperature scales of Dyck et al. ( 1996 , 1998 ) and PeiTin et 
al. ( 1998| ). CIT 3 is a regular pulsational variable with a mean 



infrared variability period of 660 days (Harvey et al. 1974 



Le Bertre 1993 Whitelock et al. 1994). The stellar surface of 
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CIT 3 is eroded by strong stellar winds with current mass-loss 
rates close to 10^^ Mo/yr. For example, various CO observa- 
tions indicate rates in the range of 0.79 up to 2.6 x 10~^ M0/yr, 
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1989 



resp. (Neri e t a. |1998| , Loup et al. |1993| , Knapp & Moms |1982 , 
Sopka et al 



1982^ . Accordingly, CIT 3 is 



Knapp et al. 

surrounded by a thick dust shell. This dust shell absorbs al- 
most all visible light radiated by the star and finally re-emits 
it in the infrared regime. The circumstellar envelope of CIT 3 
shows SiO (Dickinson et al. 119781 Desmurs et al. pOOOb, H2O 



(Dickinson 1976, Bowers & Hagen 1984) and OH maser emis- 
sion (Wilson et al. 1970, Olnon et al. 1980). The expansion ve- 
locity of the circumstellar shell amounts to 19 km/s if measured 
in the 1612 MHz OH maser lines (Baud 1198 lb, and 18 km/s to 



23 km/s, with a mean value of 21 km /s, if m easured in CO(l-O) 
and CO(2-l) transitions (Loup et al. |l993| ). 

The angular radius of the OH emitting shell region has been 
determine d by V ery Large AiTay (VLA) measurements to be 
4'.'0 (Baud |1981| ) and 4'.'4 (Bowers et al. |l983| ), resp. The cor- 
responding linear OH radius can directly be determined by the 
so-called phase-lag method utilizing the phase difference be- 
tween the red-shifted and blue-shifted peak in the 161 2Mhz OH 
spectrum during the regular IR and OH light variations. Due to 
the difference in light travel time between the emission from 
the front and the back of the shell, the variations of the se two 
masing regions are out of phase at Earth. Jewell et al. (198C) 



and van Langevelde et al. (1990) observed a phase lag of 25.5 d 
and 34 d, resp., resulting in a linear OH shell radius of i?oH = 
3.3 • lO^^cm and Rqh = 4.4 • IQi-^cm, resp. Combining the 
VLA observations with the phase-lag measurements provides 
an independent method for determining the distance to CIT 3. 
For instance. Baud ( 1981 ) derived d — 570 pc using the phase 
lag of Jewell et al. ( 1980| ), the phase lag of van Langevelde 
et al. (1990) gives d ~ 740 pc. The somewhat larger angular 
OH radius of Bowers et al. (1983 1 result in correspondingly 
smaller distances, i.e. 520 pc and 670 pc. These direct determi- 
nations can be comp ared w ith other, indirect, methods. For in- 
stance, Hyland et al. (1972) give a distance of 510pc based on 
the measured flux at maximum light (1.1 • IG^^Wm^^) and on 
the assumption that the maximum stellar luminosity amounts 
to 10**Lq. Adopt ing L = 2 • IO^Lq, would result into 720pc. 
White lock et al. (1994) applied period-luminosity relations to 
obtain the absolute bolometric magnitude which, together with 
the apparent magnitude, finally leads to d — 820 pc. In sum- 
mary, the distance to CIT 3 appears to be in the narrow range 
from 500 to 800 pc. 



Zappala et al. ( 1974 ) measured the size of CIT 3 at 2.2 /im, 
10 fim, and 20 //m by means of lunar occultations. They fitted 
the disappearence curves at 2.2 ^m and 10 fim with two-disks 
fits. At 2.2 /im, they obtained a small component with an an- 
gular diameter of 13 mas contributing 65% of the flux and a 
larger one with an diameter of 65 mas contributing 35% of the 
flux. The smaller component was interpreted to be the central 
star whereas the larger one represents the inner region of the 
dust shell. At 10 fj,m, the fit yielded one component of approx- 
imately 65 mas diameter providing 15% of the flux, and one 
of 135 mas providing 85% of the flux. In summary, CIT 3 ap- 
peared to consist of a central star with a diameter of 13 mas, 
surrounded by a dust shell with typical diameters of 65 mas at 
2.2 /im and 135 mas at 10 /im. 

Recent mid-infr ared v isibility measureme nts we re caiTied 
out by Sudol et al. ( |l999| ) and Lipman et al. ( |200C| ). Sudol et 



al. ( 1999| ) quote an upper limit of 0'.'63 for the 1 1 /im shell di- 
ameter since CIT 3 was not resolved in their 2.3 m telescope 
images. Interferometric measurements of Lipman et al. ( 200C| ) 
with the Infrared Spatial Interferometer (ISI) clearly resolve 
the object at 1 1 /im indicating an inner dust-shell diameter of 
66 mas by means of corresponding radiative transfer models. 

There are several pieces of evidence that the dust shell 
of CIT 3 deviates from spherical symmetry. Marengo et al. 
( 199S ) reported evidence for an elongation at 8.55 /im based on 
1.5 m telescope data. Neri et al. ( 1998( ) fitted their CO observa- 
tions with two-component Gaussian visibility profiles consist- 
ing of a circular and an elliptical component. The circular flux- 
dominant component coiTesponds to a spherical envelope with 
a 29'.'6 diameter, the elliptical one to an axisymmetric envelope 
structure with major and minor axis of 9'.'8 and 6'.'8, resp, and 
a position angle of —45° (E to N). 

Finally, near-infrared polarimetry (Kruszewski & Coyne 



1976, McCall & Hough 198C ) revealed that CIT 3 is highly po- 



larized in the / and J-band with polarization degrees of up to 
8-9%. The degree of polarization depends also on the phase. 
In the H- and K'-hand polarization patterns are less dominant 
(degree of polarization of the order of 1%). The high degree of 
polarization either indicates an alignment of the dust particles 
or an asymmetrical dust-shell structure. 

In this paper we present the first bispectrum speckle- 
interferometry observations of CIT 3 at 1.24/im (J-band), 
1.65 fim (H-hmd), and 2.12 /im (A"'-band). In the K'- and H- 
band, the reconstructed images are diffraction-limited, the res- 
olution in the J-band is 88% of the diffraction limit. The J-, 
H-, and A''-band resolution amounts to 48 mas, 56 mas, and 
73 mas, resp. Radiative transfer calculations have been con- 
ducted to model the spectral energy distribution as well as the 
various visibility functions in the near and mid-infrared. This 
paper is organized as follows: In Sect. ||the bispectrum speckle 
interferometric observations and data reduction are described. 
Sect. H resumes the results from available mid-infrared long- 
baseline interferometry. Sect. ^ summarizes the available pho- 
tometry, spectrophotometry and spectroscopy, and in Sect. ^ 
radiative transfer modelling and its confrontation with the ob- 
servations is discussed. Finally, a summary and concluding re- 
marks are given in Sect. ^ 

2. Observations and data reduction 

The CIT 3 speckle interferograms were obtained with the 
Russian 6 m telescope at the Special Astrophysical Observatory 
on September 19, 22 and 28, 1999. These dates correspond 
to the phases (j) = 0.74, 0.75 and 0.76 of the periodic vari- 
ability cycle of 660 days (see Sect. 4.1). The speckle inter- 
ferograms were recorded with our HAWAII speckle camera 
(HgCdTe array of 5 12^ pixels, sensitivity from I to 2.5 nm) 
through interference filters with center wavelength/bandwidth 
of 1.24/im/0.14/im (J-band), 1.65 /im/0.32 /im (iJ-band) and 
2.12 /im/0.21 fim (JC'-band). Speckle interferograms of the un- 
resolved stars HIP 7740 and HIP 6363 were taken for the com- 
pensation of the speckle interferometry transfer function. The 
observational parameters were as follows: J-band observation 
(Sept. 28): exposure time/frame 80 ms, number of frames 1963 




Fig. 1. Contour plots of the reconstructed CIT 3 images at 1 .24 /im ( J-band), 1 .65 /im (JT-band), and 2.12 /im (JsT'-band) from 
left to right. The J-, H-, and K'-bwA resolution amounts to 48 mas, 57 mas, and 73 mas, resp. The images in the K'- and ff-band 
are diffraction-limited, the resolution of the J-image is 88% of the diffraction limit . Contour levels are plotted from 1.5% to 
29.5% of peak intensity in steps of 1%. North is up and east is to the left. 
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Fig. 2. Visibility functions of CIT3 at 1.24 //m (J-band), 1.65 /im (i/-band), and 2.12 jim (if -band) from left to right. The top 
panels show the azimuthally averaged visibilities whereas the bottom panels display visibility cuts at position angles of —28° 
(long axis) and 62° (short axis), i.e. along and perpendicular to the main symmetry axis of the elongated compact structure of 
CIT 3 in the J-band (see Fig. |l]). 



(1268 of CIT 3 and 695 of HIP 6363); 7J-band observation 
(Sept. 28): exposure time/frame 20 ms, number of frames 2594 
(1723 of CIT 3 and 871 of HIP 6363); X'-band observation 
(Sept. 19H-22) : exposure time/frame 25 ms, number of frames 
6496 (3016 of CIT 3 and 3480 of HIP 7740); 1.24 /im-seeing 
(FWHM) -1'.'8; 1.65/im-seeing (FWHM) -1'.'2; 2.12 ^m- 
seeing (FWHM) --1'.'3; J-band pixel size 13.4 mas, 7J-band 
pixel size 20.1 mas, if '-band pixel size 26.4 mas; field of view 
3'.'4 (A''-band) and 5'.'1 {H-, J-band). Images of CIT 3 with 



48 mas (J-band), 56 mas (ii -band), and 73 mas (Jf'-band) res- 
olution were reconstructed from the speckle interferograms 
using the bispectrum speckle-interferometry method (Weigelt 



1977, Lohmann et al. 1983, Hofmann & Weigelt 1986). The 



H- and iiT'-band images are diffraction-limited, the resolution 
in the J-band is only slightly (12%) below the diffraction limit. 
The modulus of the object Fourier transform (visibility func- 
tion) was determined with the speckle interferometry method 



(Labeyrie 1970) 
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Figure [l| shows the contour plots of the reconstructed CIT 3 
images at 1.24 /im, 1.65 /xm, and 2.12/im. Contour levels are 
plotted from 1.5% to 29.5% of peak intensity in steps of 1%. 
The J-band image is clearly elongated along a symmetry axis 
of position angle —28°. More precisely, two structures can 
be identified: a compact elliptical core and a fainter north- 
western fan-like structure. The eccentricity of the elliptical 
core, given by the ratio of minor to major axis, is approxi- 
mately £=123 mas/154mas=0.8. The full opening angle of the 
fan amounts to approximately 40°, i.e. it extends from position 
angles —8° to —48°. In the _ff-band the elongation of the core 
is hardly visible, and the fan structure is much weaker. Finally 
in the A' '-band, the fan structure has almost disappeared be- 
ing only present on the 2.5% intensity level, and CIT 3 looks 
almost spherically symmetric. 

Figure ^ displays the coiTesponding visibility functions of 
CIT 3 at 1.24 /im, 1.65 //m, and 2.12/im, i.e. the azimuthally 
averaged visibilities as well as visibility cuts at position angles 
of —28° and 62°, i.e. along and perpendicular to the main sym- 
metry axis of the elongated compact structure and the north- 
western outflow of CIT 3 in the J-band (see Fig. |l]). Towards 
the diffraction limit, the azimuthally averaged visibility func- 
tions decline to ^ 0.2, almost reaching a plateau, at 1.24 /im, 
to ^ 0.25 at 1.65 /im and to ^ 0.45 at 2.12 /im. In principle, 
but most obvious in the J-band, all visibilities show the same 
behaviour: a sharp decline at low spatial frequencies (< 3-4 cy- 
cles/arcsec) and an extended run of the curve of smaller slope 
at larger spatial frequencies almost up to the diffraction limit. 
This indicates the existence of a compact and of an extended 
dust-shell component. 

The visibility cuts along and perpendicular to the main 
symmetry axis show similar features. Additionally, as ex- 
pected, the visiblity cuts along the symmetry axis decline 
stronger than those perpendicular to it, evidencing a larger ob- 
ject extension in this direction. The typical differences between 
the visibility cuts at large spatial frequencies amount up to 0.15 
at 1.24/im and 1.65 /im and ^ 0.1 at 2.12/im. In the J- and 
7J-band this position angle dependence can be considered as 
significant, whereas in the JC'-band the differences are of the 
order of the eiTor bar and thus are negligible. 



3. Mid-infrared visibility 

Mid-infrared visibiUties provide additional pieces of spatial in- 
formation since they are sensitive to the location of the outer, 
cooler regions of the dust shell governed by thermal dust- 
emission, whereas near-infrared visibilities are indicative for 
the inner, hot dust-formation zones and are subject to scatter- 
ing. The interferometric 11.15//m observations of Lipman et 
al. (2000) with the Infrared Spatial Interferometer at baselines 
of 4m, 9.6m and 16m are shown in Fig. || Note, that the data 
is not coeval but refers to variability phases (see Sect. 4.1) of 
0.58-0.66 (4m), 0.71-0.76 (9.6m) and 0.66-0.71 (16m). 



4.1. Variability 

CIT 3 is a regular variable with a mean infrared period of 660 d. 
Infrared lightcurves are given in Harvey et al. (1974 ), Le Bertre 
(1993), and Whitelock et al. (1994), and radio lightcurves in 
Herman & Habing (1985). The period appears to be rather 
stable with time and has not change since the 1970's. It de- 
pends slightly on the wavelength, being somewhat shorter at 
longer wavelengths. For instance, Le Bertre (1993) determined 
P = 668, 661 and 654d at A = 1.24/im, 2.19 /im, and 



4.64 /im, resp., and Herman & Habing (1985) found P = 632 d 
at A = 18 cm. The amplitudes of the brighness variations 
show, however, a st rong w avelength dependence in the near- 
infrared. Le Bertre ( |l993| ) gives Am = 2.83, 1.79, and 1.29 
at 1 .24 /im, 2.19 /im, and 4.64 /im, resp. Then, at longer wave- 
lenghts the amplitude seems to stay roughly constant. Harvey 
found Am — 1.5 at 10 /im and Herman & Habing ( 1985| ) 



Am = 1.2 at 18 cm. 

The spectrum of CIT 3 shows a 9.7 /im silicate feature in 
partial self-absorption (see Fig. H). Monnier et al. (1998) pre- 
sented multi-epoch 8-13 /im spectrophotometry (1994-1997) 
covering more than one variability cycle, and showed that the 
spectral shape of the silicate feature changes during the cy- 
cle. The brightness variation between minimum and maximum 
light was ^ 1.2 mag. 

Thus, the variability-phase dependence of the spectral en- 
ergy distribution (SED) can be characterized by wavelength in- 
dependent variations of Am ^1.2 for A ^ 10/im, a chang- 
ing spectral shape of the 9.7 /im silicate feature, and significant 
spectral shape changes for A ^ 5/im emphasizing the need of 
coeval data to construct the SED. 

All variability phases given in this paper refer to a period of 
P = 660 d and to = (or 1; maximum light) at JD 2446990 
(July 13, 1987) as in Le Bertre ([l993|). 



4.2. Data compilation 

To construct the SED various observations are available, as 
e.g . pho tometry from the near to the mid-infrared (Ulrich et 



al. 1966, Wisniewski et al. 



& Simon 1973,Dycket al 



1967, Hyland et al. 1972, Morrison 



1974, Strecker & Ney 1974, Merrill 



& Stein |1976| , Mc Carth y et al. |1977|, Epc hstein et al. 
Price & Murdock 1976, Lockwood 1985 Jones et al. 



198C 



199C 



4. Spectral energy distribution 



Fouqueet al. 1992 , Xiong et al. 1994), sub-millimeter (Sopka 
et al. 1985 ), milli meter (Walmsley et al. 1991) and radio (e.g. 
Olnon et al 1980) observations, spectrophotom etry a nd spec- 
troscopy (Lambert et al. 1980, Monnier et al. 1998, Lancon 
& Wood 2000 ), as well as photometry and spectroscopy ob- 
tained by IRAS (InfraRed Astronomical Satellite; Point Source 
Cataloque |1985| , Volk & Co hen |1989| ) and ISO (Infrare d Spac e 
Observatory; Sylvester et al. |1999[ Markwick & Millar |2000| ). 

Figure |^ displays those observations either taken close to 
or being adustable to the pulsational phase cj) of the present 
speckle measurements covering a wavelength range from 1 /im 
to 1.3 mm. This concerns the photometry of Jones et al. (1990) 
[0 = 0.77], the ISO SWS [0 0.22] and LWS [0 = 0.49] 
spectroscopy of Sylvester et al. (1999), the 400 /im data of 
Sopka et al. (1985) [(j> = 0.89], and the 1.3 mm data of 
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Walmsley et al. (1991) [<j) = 0.73] 
according to Sec 



4.1 



Phases were calculated 
Additionally IRAS fluxes (scaled to 
ISO data) are shown. Since the variability's brightness ampli- 
tude is ~ constant at mid-infrared and longer wavelenghts, the 
LWS spectrum taken at minimum light can be adjusted (fac- 
tor 1 .4, see Sylvester et al. 2000) to meet the S WS spectrum. 
Finally the complete ISO spectroscopy was scaled to match the 
ground-based photometry of Jones et al. (1990) between 9fim 
and 18/^m. The inlay of Fig. ^shows the mid-infrared spectrum 
at 4> = 0.75, interpolated from the multi-epoch spectropho- 
tometry of Monnier et al. (1998), as suitable for the present 
phase of the speckle observations, in comparison with the ISO 
data. Concerning the silicate feature, absorption is more promi- 
nent at minimum light, accordingly the shape of the ISO SWS 
spectrum (Sylvester et al. 2000) at = 0.22 fits better to the 
miminum than to the maximum light curve. 

5. Dust-shell models 



To construct models for the dust shell of CIT 3 radiative tranfer 
cal culati ons were carried out using the code DUSTY (Ivezic et 
al. 1997). The code was expanded for the calculation of syn- 
thetic visibilities as described by Gauger et al. ( 1999 1. DUSTY 
solves the radiative transfer problem in spherical symmetry uti- 
lizing the self-similarity and scaling behaviour of IR emission 
from radiatively heated dust (Ivezic & Elitzur 1997). Due to 
the observed elongation in the J-band the assumption of spher- 
ical symmetry appears to be questionable at least on smaller 
scales. However, on larger scales and for A > 2.1 /im the dust 
shell does not deviate significantly from spherical symmetry 
and its main properties can be believed to be inferable from 
such models in fair approximation. Nevertheless, it should be 
kept in mind that the derived dust-shell properties are subject 
to the geometry and are, therefore, of tentative character. 

The model construction requires the specification of the 
central radiation source, the dust properties (chemistry, grain 
sizes), the temperature at the inner dust-shell boundary and the 
relative thickness of the dust shell (ratio of outer to inner shell 
radius), the density distribution, and the total optical depth at a 
given reference wavelength. 

The models are confronted with the observed SED and the 
near- and mid-infrared visibilities. Preference is given to the 
observations close to the phase of the present speckle observa- 
tions ((j) — 0.75). Concerning the SED, we take into account 
the ground-based photometry of Jones et al. (1990) [1.6/im- 
18/im], the spectrophotometry of Monnier et al. (1998) [sil- 
icate feature], the ISO data of Sylvester et al. (1999) [5/im- 
100/xm] scaled to the ground-based photometry, as well as 
400^m (Sopka et al. |1985| ) and 1.3mm (Walmsley et al. 1991) 
data. Interstellar extinction have only minor effects (Alvarez et 
al. 2000) and was neglected. For instance, Xiong et al. (1994) 
gives for the interstellar reddening Ay = 0.18 mag referring to 
van Herk (|1965|). 



Previous (spherically symmetric) radiative transfer mod- 
els we re communicated by, e.g., Rowan-Robinson & Harris 
(|1983|), Bedijn (|l987|), Schutte & Tielens (|l989|), Le Sidaner 



the SED (from 1 /im to 1 mm) and visibilities at different wave- 
lengths into account. 

5. 1. Standard single-shell models 

An extensive model grid was calculated for standard uniform- 
outflow shells, i.e. for density distributions ^ 1/r^ within the 
dust shell. The central star was approximated by black bodies, 
however, model atmospheres (Aringer & Loidl, private comm.) 
were also taken into account. Due to its spectral type of M9- 
10, CIT3's effective temperature can be estimated to lie in 
the range of ~2500-2700K if extrapolated from t he effe ctive 
temperature scal e for giants given by Perrin et al. (1998) (see 
also Dyck et al. 1998 van Belle et al. 1999 ). However, long- 
period variables can be substantially cooler than giants of the 
same spectral type (van Belle et al. 1999), therefore models 
were calculated for Toff = 2000 to 3000 K. Concerning the 
optical dust propertie s we c onsidered differe nt set s of silicate 
grains (Drai ne & L ee 1984, Ossenkop f et al . 1992 , Le Sidaner 
& LeBertre |1993[ David & Pegourie |1995| ). In the foll owing 
sections we will mostly refer to the Ossenkopf et al. ( 1992| ) 
grains which prove to give good fits to the silicate feature. 

For the grain-size distribution we took into account grain- 
size distributions according to Mathis et al. ( 1977 , hereafter 
MRN), i.e. n(a) - fl-^-^, with 0.005 ^im < a < (0.13 
to 0.45) /im as well as single-sized grains with a = 0.05 to 
0.3 /im. Grain-size distributions were accounted for by the syn- 
thetic grain approximation which first calculates the absorption 
and scattering effiencies for each grain and then determines 
average cross sections by averaging over the size-distribution. 
The shell thickness was fixed to yout ~ '"out/^'i = 10"^, with 
Tout and ri being the outer and inner radius of the shell, re- 
spectively. Larger Fout (e.g. 10"*) affect the SED only slightly 
beyond 100/im. The model paramaters are completed by the 
dust temperature Ti at the shell's inner boundary ri , which was 
varied between 400 and 1500 K, and the optical depth, r, at a 
given reference wavelength, Arcf • Refering to Arcf = 0.55 /im, 
optical depths between 10 and 90 were considered. 

The inspection of the model grid reveals that a simultane- 
ous match of all observed properties does not appear to be pos- 
sible within the considered parameter space. In particular, it 
turned out that uniform outflow models cannot provide enough 
flux in the far-infrared and simultaneously match the observa- 
tions at shorter wavelengths. This may indicate higher densities 
in the outer shell regions than given by the l/r^ density distri- 
bution, and thus a more shallow density decline or the existence 
of regions of enhanced density due to a previous period of in- 
creas e d ma ss loss (see, e.g., Suh & Jones 1997; Blocker et al. 



& LeBertre ( |199q ), and Lipman et al. (^OOOj). The present ra- 
diative transfer models take for the first time a broad range of 



1999t |2001|). This is in line with the results of Lipman (2000) 
who found that the ISI visibility data is better fitted by a l/r^'^ 
density fall-off than by a uniform outflow. 

The best models found are shown in Fig. ^ and refer to 
an optical depth of rv=30 given by the match of the silicate 
feature. The grain sizes are strongly contrained by the near- 
infrared visibilities, and a grain-size distribution appeared to 
be better suited than single-sized grains. Considering a grain- 
size distribution, the maximum grain-size was determined to 
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be Cniax = 0.25 /im as in the original MRN distribution. The 
temperature at the inner rim of the dust shell was found to be 
in the range of 850-950 K with the best model at 900 K for the 
best match of visibilities. 

The temperature of the central star is not well contrained by 
the SED or the visibilities. For instance, black bodies of 2000 K 
and 2500 K yielded similar overall fits given the observational 
and theoretical uncertainties. However, the lunar occultation 



observations of Zappala et al. (1974), yielding a central-star 
diameter of 13 mas and a typical dust-shell diameter of 65 mas 
at 2.2 /im, seem to favour lower temperatures if one compares 
the dimensions of the respective models. For instance, adopt- 
ing Toff = 2000K abolometric flux of /boi = 1 • IQ-^Wm^^ 
determined by the SED fit gives a central-star diameter of 
8, = 13.7 mas. This leads for Ti — 900 K, to an inner dust 
shell diameter of Oi = 65.8 mas coiTesponding to ri —4.8 i?* 
in reasonable agreement with the observations. On the other 
hand, a star with Toff — 2500 K leads to a smaller central star 
of 8* = 8.8 mas and to a dust shell with 8i = 78.2 mas and 
ri = 8.9 i?*. Even though this comparion still lacks correc- 
tions due to different pulsational phases, it indicates that CIT 3 
should be cooler than, say, 2500 K. Taking into account the 
model fits and the above constraint, the best suited model is that 
with Teff = 2250 K giving 8* = 10.9 mas and 8i = 73.3 mas 
(ri = 6.7 i?*). 

The phase dependence of the silicate feature can be seen 
in the inlay of Fig. ^ showing ground-based photometry and 
spectrophotometry (Jones et al. 1990, Monnier et al. 1998) and 
the ISO data (Sylvester et al. 1999). The difference between the 
spectrophotometry (interpolated for 0=0.75) and the ISO data 
(0=0.22) becomes most obvious between 10/im and 13 /im. 
Since the scaled ISO data meets the photometry of Jones et al. 
(1990) (0=0.77) both at 9 /im and 18 /im, the phase dependence 
of the spectral shape appeared to be restricted to the above re- 
gion, at longer wavelengths flux differences are independent of 
the wavelength. The radiative transfer calculations show that 
for the considered parameters the spectrophotometry can be 
matched only in its exact shape on the expense of worsening 
the fits of other observational quantities as, for instance, the 
mid-infrared visibility. Additionally, the flux at 18/im cannot 
be matched and the flux deficit in the far infrared is consider- 
ably enlarged. 

The near infrared part of the SED tends to require some- 
what higher optical depths than, for instance, the silicate fea- 
ture or the visibilities leading to somewhat higher fluxes than 
observed. Regarding the observations this may be attributed to 
the fact that the data is for the same pulsational cycle phase but 
for different epochs, i.e. to changes of the flux variations dur- 
ing a pulsational cycle over the last decade. With respect to the 
models, the assumption of spherical symmetry has to be ques- 
tioned at shorter wavelengthts as obvious in the J-band image 
of Fig. |l]. Since the spectra of cool giants can deviate consider- 
ably from that of a black body in the optical and near-infrared, 
we used also model atmopheres (Aringer & Loidl, priv comm., 
cf. Aringer et al. 1997) instead of simple black bodies for the 
central source of radiation. The coolest atmosphere used refer 
to Toff = 2700 K and log g = 0. Due to the high optical depths, 
the SED fits remained basically unchanged, and the visibil- 



ity curves run only somewhat higher than in the coiTespond- 
ing black body calculation. For lower effective temperatures 
as, e.g., for Tcff — 2250 K being appropriate for CIT 3, these 
differences can be expected to become larger due the increas- 
ing strength of molecular absorption bands (e.g. water vapor) 
which can also be phase dependent (see, e.g., Matsuura et al. 
2001). However, due to the lack of model atmospheres in this 
regime we stayed with the black body approximation. 

We note that, while the J- and if '-visibilities decline with 
almost constant slope at higher spatial frequencies, the ii-band 
visibility shows a significant drop beyond 15 cycles/arcsec. 
This drop cannot be reproduced by any of the models. Though 
still not being significant, the development of a similar trend 
seems to be suggested by the JC'-band data close to the diffrac- 
tion limit. Such a visibility shape can be expected if, for in- 
stance, the central star is partially resolved, which, however, 
requires a stellar diameter of the order of 50 mas and certainly 
does not apply to CIT 3. The reason for this drop has to remain 
open yet but may be related to the asphericity of CIT 3. 

5.2. Different density distributions 

While the above uniform outflow models seem to match the 
observations reasonably well from the near-infrared to the mid- 
infrared , they fail to give sufficiently high fluxes beyond, say, 
20 /im. This indicates, for instance, either a flatter density dis- 
tribution than p 1 /r^ as given by the standard uniform out- 
flow model or the existence of regions with increased density 
in the outer shell due to periods of higher mass-loss in the past. 
We recalculated the main body of the grid with density distri- 
butions ranging from p ^ to p ^ and combina- 
tions of it. It turned out that distributions being generally more 
shallow can match the silicate feature and the mid-infrared vis- 
ibility but fail to fit the near-infrared visibilities. Therefore, 
two-component shells were considered with an inner uniform- 
outflow component and an outer component with a flatter dis- 
tribution, i.e. with p ^ and x < 2. The transition between 
both density distributions is at the relative radius Y2 — 7'2/ri. 
Models were calculated for Y2 = 2.5 to 30.5. 

We also calculated some uniform outflow models with den- 
sity enhancements or superwind components, i.e. with regions 
in the shell where p ^ A ■ 1 jr^ with superwind amplitudes A 
ranging from 10 to 100. However, the spectral index of the far- 
infrared tail of the SED and the visibilities seem to favour the 
above flat density distribution in the outer shell region rather 
than a superwind component. 

In order to preserve the match of SED and visibilities up to 
the mid-infrared the transition to the flatter density distribution 
should take place in sufficient distance from the inner rim of the 
dust shell. On the other hand, in order to provide sufficiently 
high fluxes in the far-infrared the exponent in the density dis- 
tribution has to sufficiently small. The best models found are 
those with p ~ 1/r^ for Yi < 20.5 and p ^ l/r^-^ for 
Y2 > 20.5. The temperature Ti at the inner rim of the dust shell 
ranges again between 850 and 950 K. Model fits are shown in 
Fig. ^. Since the transition radius Y2 — 20.5 corresponds to 
a temperature of T2 — 163 K, mostly the outer cool regions 
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Fig. 4. Fractional contributions of the emerging stellar radia- 
tion, the scattered radiation, and of the dust emission to the to- 
tal flux as a function of the wavelength for a two-component 
model. The density decreases ^ Xjr^ for Y2 < 20.5, and 
~ I/t-I-S fQj- > 20.5. (see Fig. ||). Model parameters are: 
black body. Toff = 2250 K, Ti = 900 K, rcss^m = 30, 
Ossenkopf et al. (1992) silicates, and a Mathis et al. (1977) 
grain size distribution with flmax = 0.25 /im. 



of the dust shell are concerned and the spatial dimension of 
the dust shell changed only slightly. For Toff — 2250 K and 
Ti = 900 K, the inner rim of the dust sheU is at ri/i?* = 6.62 
corresponding to an angular diameter of Oi ~ 71.9 mas. 

These results can be compared with the radiative transfer 
model of Lipman et al. (2000). Based on fits of the ISI 1 1 /im 
visibility data (and the silicate feature) they also found a 1 /r^-^ 
density fall-off to be better suited to match the observations. 
However, an overall l/r^^ density law appears to be in conflict 
with the present near-infrared visibilities leading to too high 
densities in the inner parts of the shell. The present modelling 
suggests, as mentioned above, the presence of a sufficiently 
broad uniform outflow area in the inner shell region before the 
flatter density distribution takes over. 

Fig. ^ shows the fractional contributions of the emerging 
stellar radiation, the scattered radiation, and of the dust emis- 
sion to the total flux as a function of the wavelength for the 
above two-component model. Whereas the radiation at 1 1 /im 
is almost completely given by thermal dust emission, the radia- 
tion in the near-infrared domain is strongly determined by scat- 
tering. The contribution of scattered light amounts to 20.6%, 
42.3% and 79.2% in the if'-, H- and J-band, resp., while 
thermal dust emission accounts for 21.2%, 4.8% and 0.2%, 
resp., and the emerging direct stellar light to 58.2%, 52.9%, 
and 20.6%. resp. It demonstrates that for the given grain-size 
distribution 49.3%, 89.8% and 99.7% of the dust-shell emis- 
sion is due to scattering. We note that these derived attributes 
of the near-infrared emission are subject to some uncertainties 
due to their dependence on the dust properties and the assump- 
tion of spherical symmetry. For instance, though the grain-sizes 
are strongly constrained by the near-infrared visibilities, prin- 
cipal ambiguities remain due to, e.g., the used optical dust con- 
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Fig. 5. Normalized intensity at 1.24 /im, 1.65 /im, 2.12 //m, and 
11/im (from top to bottom) vs. angular displacement -d for a 
two component model with a uniform outflow in the inner dust 
shefl region and p ~ 1/r^-^ for Y2 > 20.5 (see Figs. ||-||). 
The central peak belongs to the central star with a radius of 
5.5 mas. The inner hot rim of the circumstellar shell has a radius 
of 35.9 mas, the cool far-out component with the flatter density 
distribution is located at a radius of 737 mas beyond the scale 
of this figure. In the if' -band the inner rim of the circumstellar 
shell is noticeable limb-brightened. The model refer to Toff = 
2250 K, Ti = 900 K, and T(0.55/im) = 30. The calculations 
are based on a black body, Ossenkopf et al. (1992) silicates, 
and a Mathis et al. (1977) grain size distribution with Omax — 
0.25/im. 



stants depending on shape and chemistry of the grains. Since 
the resolution of the present speckle observations is of the or- 
der of the inner dust-shell boundary, future near-infrared in- 
terferometric measurements of higher resolution can yield fur- 
ther constraints on the models and dust properties, resp., by al- 
lowing a direct measurement of, e.g., the distribution of direct 
stellar and scattered radiation. The corresponding spatial inten- 
sity distributions at 1.24 /tm, 1.65 fj,m, 2.12 /im, and 11/im are 
shown in Fig. |[ The optical depths are T(1.24/im) = 5.1, 
T(1.65/im) = 3.1, r(2.12/im) = 1.9, andr(ll/im) = 2.2. 
Although exposing an optical depth larger than unity, the inner 
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rim of the circumstellar shell is noticeably limb-brightened in 
the A''-band. 

The luminosity of CIT 3 is 7791 L0 if d = 500 pc and 
19945 Lq if d = 800 pc. The former value would predict CIT 3 
to be an AGB star of lower mass whereas the latter indicates 
a more massive AGB star with a core mass larger than, say, 
0.8 Mq. Adopting a dust-to-gas ratio of 0.005 and a specific 
dust density of 3gcm^'', the present-day mass-loss rate can 
be determined to be M = 1.3 • 10~^Mo/yr if d = 500pc 
and M = 2.1 • lO-'^Mg/yr if d ^ 800pc which is in line 
with the CO observations (see Sect. [l|). The smaller exponent 
of 1.5 in the density distribution beyond ^2=20. 5 means that 
the mass-loss rate has decreased with time in the past provided 
the outflow velocity kept constant. Adopting i;cxp=20 km/s, the 
kinematical age of the inner uniform outflow zone can be de- 
termined to be 87 yr. 

6. Summary and conclusions 

We presented the first near-infrared bispectrum speckle inter- 
ferometric observations of the oxygen-rich, long-period vari- 
able AGB star CIT 3 in the J-, H-, and K'-hmd. The observa- 
tions were taken with the Russian SAO 6m telescope. The reso- 
lution in the J-band is 88% of the diffraction-limit (48 mas), in 
the H-, and i^'-band the resolution is fully diffraction-limited 
(56 mas, and 73 mas, resp.). 

While CIT 3 appears almost spherically symmetric in the 
H- and X'-band it is clearly elongated in the J-band along 
a symmetry axis of position angle —28°. Two structures can 
be identified: a compact elliptical core and a fainter north- 
western fan-like structure. The eccentricity of the elliptical 
core, given by the ratio of minor to major axis, is approxi- 
mately £=123 mas/154mas=0.8. The full opening angle of the 
fan amounts to approximately 40°, i.e. it extends from position 
angle -8° to -48°. 

Extensive radiative transfer calculations have been caiTied 
out and confronted with the observations taking into account 
the spectral energy distribution ranging from 1 fim to 1 mm, 
our near-infrared visibility functions at 1 .24 ^m, 1 .65 fj,m and 
2.12 /im, as well as mid-infrared interferometric measurements 
at 11 /im. It has to be noted that these calculations rely on the 
assumption of spherical symmetry. Given the observed elonga- 
tion in the J-band, this assumption appears to be questionable 
at least on smaller scales. However, the main dust properties 
can be believed to be inferable from such models in fair ap- 
proximation since the deviations from spherical symmetry ap- 
pear only on smaller scales and for A < 2.1 /im. Nevertheless, 
the derived dust-shell properties are subject to the geometry 
and are, therefore, of tentative character. It is due to future two- 
dimensional radiative transfer calculations to prove the degree 
of reliance of the presented models. 

The best model found to match all observations refers to 
a cool central star with Tcft = 2250 K which is surrounded 
by an optically thick (r(0.55/im) = 30) dust shell. The low 
temperature of the central star is constrained by the models 
but in particular a lso by the lunar occultation observations 
of Zappala et al. ( 1974 ) who derived a stellar diameter of 
13 mas. For the given bolometric flux, this diameter can only 



be achieved for low effective temperatures. The present mod- 
els give a central-star diameter of = 10.9 mas being in line 
with these observations. Extrapolating the temperature scale of 
Perrin et al. ( 1998 1 based on short-period variables would pre- 
dict Toff ~2500-2700K for M9-10 stars. The lower effective 
temperature found for CIT 3 is consistent with the finding that 
long-period variables can be substantially cooler than giants of 
the same spectral type but of shorter period (van Belle et al. 



1999) 



The temperature at the inner rim of the dust shell amounts 
to Ti — 900 K. This is somewhat cooler than the usually as- 
sumed condensation temperatures for silicate grains of 1000- 
1200 K, but still well within the physical range of dust con- 
densation temperatures and in reasonable agreement with the 
results of previous radiative transfer models (as e.g., Ti 
850 K, Le Sidaner & LeBertre |1996| ; - 900 K , Bedijn [1987 
or ~ 1000 K, Rowan-Robinson & Harris 1983 Lipman et al 



200c ). The inner dust shell boundary is located at ri = 6.6i?* 
corresponding to an angular diameter of 81 = 71.9 mas. This 
angular diameter is in accordance with the 2.2 jim diameter of 
65 mas derived by lunar occultation observations (Zappala et 
al. 1974) tracing the hot dust and therefore the very inner parts 



of the dust shell. 

Grains with various optical properties were used, the sil- 
icate grains of Ossenkopf et al. (1992) turned out to be su- 
perior to the other grain types, for instance in fitting the sil- 
icate feature. The grain sizes are strongly constrained by the 
near-infrared visibilities, and a grain-size distribution proved 
to be better suited than single-sized grains. We considered 
grain-size distribution according to Mathis et al. (1977), i.e. 
n{a) ^ a^^-^, with 0.005 /im < a < Omax and determined 
ftmax = 0.25 /im as in the original MRN distribution. 

Uniform outflow models, i.e. density distributions with p ^ 
1/r^, turned out to match the near- to mid-infrared part of the 
SED and the corresponding visibilities but to considerably un- 
derestimate the flux beyond 20 /im. This indicates, for instance, 
either a more shallow density distribution than p ~ 1/r^ as 
given by the standard uniform outflow model or the existence 
of regions with increased density in the outer shell due to pe- 
riods of higher mass-loss in the past. The model grid was re- 
calculated for a large variety of accordingly changed density 
distributions, and a two-component model existing of an inner 
uniform-outflow shell region (p l/r'^) and an outer region 
where the density declines as p ^ 1/r^^ proved to be best 
suited to match the observations. A p ^ density dis- 

tribution has also been proposed by Lipman et al. (2000) by 
fits of the 1 1 /im visibility. The present modelling suggests, 
however, that the flatter density fall-off should be restricted 
to the outer shell regions which is, for instance, constrained 
by the near-infrared visibilities. The transition between both 
density distributions takes place at r2 = 20.57-1 — 135. 7i?* 
where the dust-shell temperature has dropped to T2 = 163 K. 
Since these changes concern mostly the very cool regions of 
the shell, the other so-far discussed dust-shell properties remain 
basically unchanged. Provided the outflow velocity kept con- 
stant, the more shallow density distribution in the outer shell 
indicates that mass-loss had decreased with time in the past 
of CIT 3. Adopting v^xp — 20 km/s, the termination of that 
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mass-loss decrease and the begin of the uniform-outflow phase 
took place 87 yr ago. The present-day mass-loss rate can be 
determined to be M = 1.3 • lO^^Mo/yr ii d = 500pc and 
M — 2.1 ■ 10^^ M0/yr if d = 800 pc being in agreement with 
the CO observations, and the corresponding stellar luminosity 
is 7791 L0 and 19945 L©, resp. 

CIT 3 proved to be among the most interesting far-evolved 
AGB stars due to its infrared properties. Moreover, the aspher- 
ical appearence of its dust shell in the J-band puts it in one 
line with the few AGB stars known to expose near-infrared as- 
phericities in their dust shells. The development of such as- 
phericities close to the central star suggests that it is in the very 
end of the AGB evolution or even in transition to the proto- 
planetary nebula phase where most objects are observed in 
axisymmetric geometry (Olofsson 1996). However, in contrast 
to other objects, CIT 3 shows these deviations from spherical 
symmetry still only in the J-band, which is almost completely 
dominated by scattered light. This suggests that CIT 3 had just 
started to form aspherical structures and is in this regard still in 
the beginning of its very final AGB phase. If so, CIT 3 is one 
of the earliest representatives of this dust-shell transformation 
phase known so far. 
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Fig.3. Model SED and visibilities at 1.24/xm, 1.65 ^im, 2.12/xm, and 11 fim for models with (Toff, Ti)=(2000K,850K) and 
(2500K,950K), resp. Both uniform outflow models (p ^ solid and long-dahed lines) and models with a flatter density 

distribution in the outer shell region (p ~ for Y < 20.5 and p ~ l/r^-^ for Y > 20.5; short-dashed and dotted lines) 
are shown. The optical depth ro.ss^m is 30. The calculations are based on a black body, Ossenkopf et al. (1992) silicates, and 
an MRN grain size distribution with Omax = 0.25 ^m. The symbols refer to the observations. SED — O; Jones et al. (1990) 
[(/> = 0.77]; □: Sopka et al. (1985) [(f) = 0.89]; O: Walmsley et al. (1991) [0 = 0.73]; y: IRAS fluxes (scaled to ISO data); *: 
spectrophotometry of Monnier et al. (1998) [interpolated for (p = 0.75]; thick line: ISO spectroscopy of Sylvester et al. (1999) 
[adjusted to match the Jones et al. (1990) photometry at 9/im and 18/im]. The inlay shows the Jones et al. (1990) and Monnier 
et al. (1998) data together with the ISO spectroscopy of Sylvester et al. (1999). Near-infrared visibilities — A: present speckle 
observations. Mid-infrared visibility — A: ISI observations at 4m, 9.6m and 16m baseline of Lipman et al. ( 2000| ). 



